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We apply to the calculation of the pressure of a hot scalar field theory a method 
that has been recently developed to solve the Non-Perturbative Renormalization 
Group. This method yields an accurate determination of the momentum dependence 
of n-point functions over the entire momentum range, from the low momentum, pos- 
£■"»«. ■ sibly critical, region up to the perturbative, high momentum region. It has therefore 

o. 

the potential to account well for the contributions of modes of all wavelengths to 
the thermodynamical functions, as well as for the effects of the mixing of quasi- 



particles with multi-particle states. We compare the thermodynamical functions 
obtained with this method to those of the so-called Local Potential Approximation, 
and we find extremely small corrections. This result points to the robustness of the 
quasiparticle picture in this system. It also demonstrates the stability of the overall 
approximation scheme, and this up to the largest values of the coupling constant 
that can be used in a scalar theory in 3+1 dimensions. This is in sharp contrast to 
perturbation theory which shows no sign of convergence, up to the highest orders 
that have been recently calculated. 

PACS numbers: 05.10.Cc, ll.10.Wx, 11.15.Tk 



I. INTRODUCTION 

There has been much effort devoted in the recent years to the development of finite 
temperature field theory, in particular in the context of Quantum Chromodynamics (QCD) 
whose thermal properties are central to the understanding of heavy ion collisions at high 
energy. Equilibrium properties of hot QCD are calculable from lattice gauge theory (for a 
review, see e.g. [1]), but there is a need to develop semi-analytical tools to understand the 
results of such calculations, with the hope that such tools may also allow one to approach 
non-equilibrium situations. 

Weak coupling expansions are among such tools. In the case of QCD, the use of per- 
turbation theory is motivated by the asymptotic freedom that leads to a small effective 
coupling at high temperature. However, strict perturbation theory does not work at finite 
temperature: it exhibits indeed very poor convergence properties, even in a range of values 
of the coupling where good results are obtained at T = 0. This difference of behavior of 
perturbation theory at zero and finite temperature can be understood from the fact that, at 
finite temper, tfie e X pan S1 on parameter invo.ves both the coupling and the magnitude 
of thermal fluctuations (for a recent review, see [2|; see also [3[). In that respect, the problem 
is not specific to QCD: Similar poor convergence behavior appears also in the simpler scalar 
field theory |4(, and has also been observed in the case of large- iV 4 theory [5j. Similar 
observations can be made also in the case of QED |6|]. 

In this paper we restrict ourselves to the case of the simple theory of a scalar field ip, with 
a quartic self-coupling g 2 (p^, for which high order calculations were recently completed. Thus 
the pressure is known to order g 8 \ng [7|, and the screening mass to order g 6 \ng [8[. Such 
calculations were made by exploiting effective theory techniques, in particular dimensional 
reduction, that rely on a separation of scales of various degrees of freedom in the hot scalar 
plasma: hard modes with momenta k ~ T that contribute dominantly to the pressure and 
are weakly coupled, and soft modes with momenta k ~ gT that are more strongly coupled. 
In QCD, another scale concerns the ultrasoft modes, with momenta k ~ g 2 T that remain 
strongly coupled for any value of the coupling constant. This scale, relevant only in the 
case where massless modes exist at finite temperature, does not play any role in the scalar 
field theory. The separation of scales that allows the organization of the calculation using 
effective field theory disappears when the coupling is not too small: then the various degrees 



of freedom mix and the situation requires a different type of analysis. The purpose of this 



paper is to apply to this problem t 



re Non-Perturbative Renormalization Group (NPRG) 



(for reviews of this method see e.g. SHUj) 



To do so, we shall rely on an elaborate approximation scheme that has been developed 
recently in order to obtain a good determination of the momentum dependence of the n- 
point functions [12]. This new approximation has been tested on the O(N) models, for which 
it provides excellent critical exponents, and more generally, an excellent description of the 
momentum dependence of the 2-point function, from the low momenta of the critical region, 



all the way up to the large momenta of the perturbative regime [121]. One may then expect 
this method to capture accurately the contributions to the thermodynamical functions of 
thermal fluctuations from various momentum ranges, and hence handle properly the mixing 
between degrees of freedom that takes place as the coupling grows. Since it involves also 
non trivial momentum dependent self-energies, the method also encompasses effects related 
to the damping of quasiparticles, or their coupling to complex multi-particle configurations. 
The present paper may be viewed, in its spirit and goals, as a follow up of the analysis 
presented in Ref. [14 ]. However it departs from it in two ways. The first difference is 
of a technical nature: the new approximation scheme is better justified when one uses 
an Euclidean symmetric four dimensional regulator that cuts off the contribution of high 



Matsubara frequencies. In contrast, in Ref. 14J, we used a three dimensional regulator, and 



performed analytically the (untruncated) sums over the Matsubara frequencies. The second 
difference is that we use a new, much more accurate, approximation scheme to solve the 
NPRG equations, as mentioned above: The Local Potential Approximation (LPA) used in 



Ref. [l4j can be viewed as the zeroth order in this new approximation scheme. In physical 
terms, the LPA corresponds to an approximation where the degrees of freedom of the hot 
scalar plasma are massive quasiparticles. The new scheme goes beyond that simple picture. 
As it turns out, the results obtained are not too different from those of the LPA. This 
stability of the results against improvements in the approximation suggests that the scheme 
that we are using to solve the NPRG equations may give already, at the level at which it 
is implemented here, an accurate representation of the exact pressure, and this over a wide 
range of coupling constants. It also indicates that for such a system the quasiparticle picture 
is presumably robust. 

The outline of this paper is as follows. In Sec. [Ill we present a brief introduction to 



the NPRG and the approximation scheme of Ref . [12] , indicating specific features of finite 
temperature calculations. In Sec. lIHI we integrate the flow equations numerically and discuss 
the results obtained. The last section summarizes the conclusions. In the Appendices we 
give details about the numerical integration needed to calculate the flow of the pressure, 
and we discuss specific features of the exponential regulator used in our calculations at 
finite temperature. 

II. THE NPRG IN THE BMW APPROXIMATION SCHEME 

We consider a scalar field theory with the classical (Euclidean) action 

-l/T 



S = ^ T dr|d 3 x{i(^(x)) 2 + ^^) + ^ 4 (x)} 



(2.1) 



where T is the temperature. Our goal is to calculate the thermodynamical pressure for this 
scalar theory, using the NPRG. We follow here Ref. [lOj, and add to the original action S a 
regulator 



~J R K {q)<p{qM-q), (2.2) 



2„, 

where the parameter k runs continuously from the microscopic scale A down to 0. We use 
the notation 



"i 
^2 _ , ,2 i „2 



(2tt)' 



with q = (u n , q), q — u;„ + q , and u n = 2nnT are the Matsubara frequencies. The role of 
AS K is to suppress the fluctuations with momenta \q\ < k, while leaving unaffacted those 
with |g| > k (\q\ = yep 1 )- This is achieved with a cut-off function R K (q) that has the 
following properties: R K {q) — > k 2 when \q\ <^C k, and R K {q) — > when \q\ 3> K. The precise 
form of the function R K (q) used in our calculation will be discussed later. 

The effective action T K [<f)] associated to the action S + AS K obeys the exact flow equation 



151 ] (with d t = nd K ): 

dtTM = \f d t R K (q) G K [q, -q; 0], (2.4) 

where G K [q, —q; <f)} is the full propagator in the presence of the background field 0: 

G-% -q; 0] = T^[q, -q- 0] + R K (q), (2.5) 



with T K [q, — q; 0] the second functional derivative of T K [(j)] w.r.t. <\>. The initial condition 
of the flow is specified at the microscopic scale k = A: at this point, we assume that the 
fluctuations are completely frozen by AS K , so that T k= a[0] w S[4>]. The effective action T[<f>] 
of the scalar field theory is obtained as the solution of Eq. (12.4ft for k — > 0, at which point 
R K (q) vanishes. 

When <p is constant, the functional T K [<p] reduces, to within a volume factor, to the 
effective potential V K (4>). The flow equation for V K follows from that of the effective action 
T K , Eq. (j2.4p . when restricted to a constant (p. It reads 

dtV K (p) = \fd t R K {q)GMp), (2-6) 

where 

G-\q,p) = lf(q,p) + R K (q), p = £. (2.7) 

We used here the simplified notation T ( K (q, p) in place of 1^ \q, — q, p) for the 2— point 
function in a constant background field, and similarly for G(q, p). Also, we have set p = <ft 2 /2, 
a notation to be used throughout (when <ft is constant). The pressure P is related to the 
effective potential by 

P K (T) = - [V K (T, p = 0) - V K (T = 0, p = 0)] . (2.8) 

The equation for the effective potential may be viewed as the equation for the "zero- 
point" function in a constant background field. By taking two derivatives with respect to <p 
and letting be constant, one obtains the equation for the 2-point function in a constant 
background field: 



dtT^\p) = / d t R K {q)Gl{q) 

Jq 

| r l 3) (P, Q, ~P ~ q)G K (q + p)T® {-p, p + q, -q) - -i^ 4) (p, -p, q, -q) \ 



Jq 

X 

(2.0) 



In this equation, all the n-point functions depend on the constant background field (p. This 
has not been indicated explicitly in order to alleviate the notation. 

The flow equations (12. 6ft and ( 12.91) are the first equations of an infinite tower of coupled 
equations for the n-point functions, whose solution requires some truncation. We shall use 
the truncation scheme proposed recently by Blaizot, Mendez-Galain and Wschebor (BMW) 
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12| . implemented here at its lowest non-trivial order: in this case we need only consider 
the equation for the effective potential and that for the 2-point function, that is, Eqs. (12. 6p 
and (12.91) . on which we shall perform approximations that are described next. For a more 
complete discussion we refer to Ref. [161 ]. 

The lowest level of the approximation scheme corresponds to a widely used approximation, 
usually referred to as the Local Potential Approximation (LPA). It consists in assuming that 



for all values of k the effective action takes the form 10] 



T M = J Q /d 3 xj^(90) 2 + K(0)}, (2.10) 
which is tantamount to assume that the 2-point function is of the form 

T^(q,p) = q 2 + ml(p), m 2 K (p) = %V K . (2.11) 
With this Ansatz, Eq. (J2.6P for V K becomes a closed equation that we write as follows 

d t V K (p) = ±I U (2.12) 
where I\ is one of the following integrals 

J n (p) = f d t R K {q) G K (p + q)G n K -\q), I n = J n (p = 0). (2.13) 

Jq 

Note that Eq. ( I2.12p is formally an exact equation, and it would yield the exact effective 
potential if ii were calculated with the exact propagator. In the LPA, the equation keeps 
the same form, but the propagator is given by Eq. ( 12. lip , where rn 2 .{p) is itself determined 
by the potential (thereby making Eq. ( 12. 12ft a closed, self-consistent equation). 

The next order of the approximation scheme, to which we refer to as the leading order 



of the BMW method 12J , or here simply as BMW for briefness, consists in neglecting the 
loop momentum q in the 3 and 4-point functions in the right hand side of Eq. ( 12. 9p . Once 
this approximation is made, the corresponding n-point functions can be obtained as the 
derivatives of the 2-point function with respect to the constant background field: 

T^\p, -p, 0, 0) = fylf >(p, 0), T^(p, -p, 0, 0, 0) = dl^ 2 \p, 0). (2.14) 

The equation for the 2-point function becomes then a closed equation 

d t Vf{p,p) = J a (p) (d^ 2 \p,p)) 2 - \l 2 d 2 TW(p,p). (2.15) 



There is however a subtlety: this equation for the 2-point function is coupled to that of 
the effective potential, because r« (p = 0,p) = d^V K (4>). In order to properly implement 
this coupling, we treat separately the zero momentum (p = 0) and the non-zero momentum 
(p ^ 0) sectors, and define 

r?W) = P 2 + A K (p,p) + ml(p), (2.16) 

where m 2 K (p) = d 2 V K (p) is obtained by solving the equation for the effective potential. The 
equation for A K (p, p) can be easily deduced from that for T^ 2 \ i.e., from Eq. ( 12.151) by 
subtracting the corresponding equation that holds for p — 0. It reads 

d t A K (p, p) = 2pJ 3 (p, k, p) [u K (p) + A' K (p, p)f - 2pl 3 (n, p) u 2 K (p) 

- l -h{K,p) [A' K (p,p)+2pA'>,p)], (2.17) 

where the symbol ' denotes the derivative with respect to p, and we have set u K (p) = 
dm 2 K (p)/dp. 

This equation ( 12.171) . together with Eq. ( 12.6f) for the effective potential and that for the 
propagator 

G~\q, p)=q 2 + A K (q, p) + m 2 K (p) + R K (q), (2.18) 

constitute a closed system of equations for A K (p, p) and V K ((p). This can be solved with the 
initial condition T A (p; p) = p 2 + m 2 + up, where m 2 and u are essentially (to within small 
ultraviolet cut-off corrections) the parameters of the action ( 12.11) . 

We now specify the regulator R K (q) that we have used in our calculation. We take it of 
the generic (Euclidean symmetric) form 

R K (q) = Z K K 2 r(q) } q=l, (2.19) 

K 

where Z K is a function of k only, to be specified shortly, and the function r(q) is a smooth 
function of its argument. In most of our calculations, we have used an exponential regulator 
of the form 

r(q) = -^1, q = ~, (2.20) 

e q — 1 k 

where a is a free parameter. We have also considered an alternative regulator, given by 

r (q) = ae~^ 2 -^\ (2.21) 
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with parameters (3 = 1/2 and 7 = 1/24 chosen so that the Taylor expansion around q = 
agrees with the regulator (12.201) through order 0(q A ), and leave only the prefactor a as a 
free parameter 1 . The regulators (12.201) and (12.211) are suited for the BMW approximation 
whose justification relies on both momenta and frequencies being cut-off. For comparison, 
we have also solved the LPA with these exponential regulators. We shall also compare to the 
LPA results of Ref. [14j where the Litim regulator r(q) = (1 — q 2 )6(l — q 2 ) [l7j| was used in 
integrals over three momenta, the Matsubara frequencies being integrated over analytically. 

In the absence of any approximation, the physical quantities such as the pressure should 
be strictly independent of the cut-off function, and in particular of the value of the parameter 
a that we have introduced in the regulators (I2.20p and ( I2.2ip . In practice, we find a (weak) 
residual dependence on the parameter a, and a study of this spurious dependence provides 
an indication of the quality of the approximation. It should be emphasized however that 
this dependence on a is rather small, as we shall see, and we have made no effort to perform 
a systematic study of the dependence of our results on its value. Nor did we explore the 
effects of enlarging the space of possible variations by allowing for instance the parameters 
/3 and 7 in ( 12.211) to take arbitrary values. 

The factor Z K in Eq. (I2.19p reflects the finite change in normalization of the field between 
the scale A and the scale k. It is defined by 

z - d 

op 

where po is a priori arbitrary but chosen here to correspond to the minimum of the effective 
potential. This factor Z K enters also the definition of the dimensionless variables that are 
used in the numerical solution. Thus, for instance we define 

q = nq, p = Kn 2 Z~ l p, m 2 K = Z K K 2 m 2 K , u K = Z 2 K K~ x u K , (2.23) 

with m K = rrik(p — 0), u R = u K (p), and K = l/(327r 2 ). In the next section, when we discuss 
the results obtained by solving numerically Eqs. (12.171) and (I2.6p . we shall set g 2 = u K /24. 



o-4 2 W) 



(2.22) 

P=0,p=PQ 



1 We are grateful to B. Delamotte and H. Chate for suggesting the usefulness of this regulator 



III. RESULTS 

To calculate temperature dependent physical quantities, we need to evaluate the flow of 
these quantities at both zero temperature and at finite temperature. We follow here the 

nn 

strategy exposed in Refs. [14], H8j . Conceptually, one starts with given physical parameters 
at zero temperature at the scale k = 0. One then removes quantum fluctuations step by 
step by integrating the flow equations upwards from k = to A in order to arrive at "bare 
quantities" at a chosen scale A. If A is chosen big enough, only renormalizable local operators 
survive in the effective action, and the system can then be described by a simple set of bare 
parameters g\, m\, and Z\, Starting from these bare parameters one then follows the flow 
downwards from k = A to 0, but this time with the temperature T turned on. The physical 
quantities are then obtained at K — 0. 

In practice, for reasons of numerical stability, the flow is never integrated upward, but 
always downward, for both zero and finite temperature. For a given bare coupling g^ (and 
Z\ = 1), one adjusts the bare mass m\ at the scale k = A by bisection, so as to arrive at 
a vanishing mass itlq = at the end (n = 0) of the zero temperature flow, to the desired 
accuracy. Keeping the same bare parameters, one then turns on the temperature and runs 
the flow again. Note that we consider specifically here a massless theory in order to be able 
to compare our results with the known results of high order perturbation theory. But the 
same analysis could be done for any value of tuq. 

Since the zero-temperature coupling constant vanishes at k = for any g^ (the so-called 
triviality of (p A theory in d — 4), one adjusts the coupling constant at a finite scale k > 0. 
In accordance with what is commonly done in perturbation theory, we choose to do that at 
k = 2ttT on the T = flow. This procedure induces a specific scheme dependence, attached 
to the choice of the regulator, which should be taken into consideration when comparing with 
results that are obtained in another scheme, for example the minimal subtraction scheme in 
perturbation theory. 

In order to keep the zero-temperature flow under control, we introduce dimensionless 
variables [10(] . Although these are not optimal for the thermal flow, which freezes in dimen- 
sionful variables (but diverges in dimensionless variables) as n — > 0, it is still advantageous 
to use the same scaling of variables as for the zero-temperature flow in order to achieve 
high-precision cancellations between zero and finite temperature flows for larger /t>T. 
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Figure 1: (Color online) Flow of various quantities in the BMW approximation as a function 
of the flow parameter n at zero and finite temperature. The left column shows results at small 
bare coupling (g\ = 0.88) while the right column shows results at large bare coupling (g\ = 2.8). 
Depicted are the flows of the quartic coupling, the mass, the pressure (Po = 7r 2 T 4 /90), and the 
Z-factor (top to bottom). Vertical lines indicate the positions where n = gT, T, and 2-kT. 



We have adapted to finite temperature the numerical strategy that has been deve 
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oped 



to solve the BMW equation at zero temperature in the context of O(N) models [13|, [16]. 
This method puts propagator and potential on a regular two-dimensional grid in |q| and p 
variables and solves the flow equations using the Euler method. At finite temperature, the 
flow equation (12.61) involves a summation over Matsubara frequencies. This summation is 
infeasible in practice at large values of k, but at sufficiently large k>T, the corresponding 
vacuum integral constitutes the dominant contribution and any thermal contribution is 
exponentially suppressed. Practically, a switching temperature T switch ^> T is introduced 
such that, at k = T switc h, the numerical code switches from vacuum integration to thermal 
integration with Matsubara frequencies. At this point, the values of the vacuum propagator 
G K (q,p) which are stored on a two-dimensional grid (\q\,p) have to be interpolated in order 
to yield values on a three-dimensional grid (u, |q|,p) using cubic splines. Good results have 
been achieved on a n q x n p = 90 x 50 grid with g max = 9 and p = 5 (in dimensionless variables, 
see Eq. (I2.23[) ) with 30 Matsubara terms at the switching temperature T switc h which is varied 
in the range AtcT < T sw i tc h < QttT. To speed up the calculation, the number of Matsubara 
terms can be gradually reduced as the flow proceeds, without affecting accuracy. Further 
details on the integration of the pressure flow are given in Appendix A. 

Figure [1] shows the flow of various quantities (coupling constant, mass, pressure and Z- 
factor) in the BMW approximation. (The pressure Po is that of the non-interacting scalar 
plasma, i.e., Po = tt 2 T 4 /90.) These flows are remarkably similar to the corresponding ones 
obtained within the LPA approximation in Ref. [14J. At finite temperature, the flow starts 
to deviate from the vacuum flow between 2nT and T, and stabilizes shortly below the scale 
gT (with g = g(2irT)). Note that gT is the leading order value at small g of the thermal 
mass of the excitations (see e.g. Fig. [2]): when k reaches this value the mass takes over the 
role of an infrared cut-off, which freezes the flow. For the bare value g\ = 0.88 the flow of 
the coupling reaches a value g(n = 2tcT) = 0.66, while for the larger bare coupling g\ = 2.8 
the value is g{2irT) = 1.1. Vertical lines indicate the positions of k = gT, T, and 27rT for 
the temperature T = 10 _2 A. The bottom figures displays the flow of the Z-factor, which 
is trivial in the LPA approximation (where Z = 1 for all values of k). It turns out that 
this factor deviates only moderately from 1 in the BMW approximation. This indicates 
that the mixing of single particle excitations with multi-particle states is very mild, at least 
within the BMW approximation, and suggests that the quasiparticle picture is robust in this 
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Figure 2: (Color online) Mass and pressure as function of the coupling. The various resummation 
and RG methods compared are 2 PI, LPA with exponential regulator (exp), LPA with Litim 
regulator, and BMW. Shown are also perturbative results through order g^logg for the mass and 
g 8 log g for the pressure. The g 7 and g 8 log g curves for the pressure almost lie on top of each other, 
as do the BMW and LPA (exp) curves for the mass. 

system. This could be confirmed by a complete calculation of the single particle spectral 
function, which is beyond the scope of this paper. 

Figure |2] shows a comparison of different approximations to the mass and the pressure 
as a function of the coupling g(2nT). To obtain these plots, results similar to the ones 
shown in Fig. [H and obtained for various bare couplings and ratios T/A (where values 
up to T/A = 1/10 have been used in order to extend the plot range to large values of 
g(2nT)), have been combined into single plots. Results obtained within the LPA and the 
BMW approximations to the NPRG are compared to perturbation theory through order 
g 6 \ogg for the mass |8j and q s log q for the pressure [7|, and also to the result of the 2- loop 



2PI resummation from Ref. 14|. The BMW results were all obtained with the exponential 
regulator (12.20p . while the LPA results were obtained with this same exponential regulator, 
as well as with the 3-dimensional Litim regulator used in Ref. |14j. For a given choice of 
regulator, one sees in Fig. [2] that the difference between the LPA and the BMW results for 
both the mass and the pressure is tiny: one does not gain much, for the thermodynamics, in 
improving the treatment of long wavelength degrees of freedom and incorporating explicit 
momentum dependence in the self-energy The stability of the BMW approximation scheme, 
comparable to that of the 2PI resummation, should be contrasted with the wild oscillatory 
behavior of the successive orders of perturbation theory. 
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The plots in Fig. [2] depend on the choice of the prescription for obtaining the coupling at 
a particular scale (in this case at the scale \i = 2ttT), and hence on the regulator: depending 
on the regulator, a given value of the coupling constant g(2irT) will correspond to different 
bare coupling constants g&, and hence to different initial conditions for the temperature flow. 
Together with the obvious dependence of the flow itself on the regulator this will affect the 
final result. The regulator dependence that is visible in the LPA results in Fig. [2] is to be 
attributed in part to the way the results are plotted, namely as a function of g(2irT). Most 
of such a "scheme dependence" can be eliminated by plotting only physical quantities. This 
is realized in Fig. [3] which shows the pressure as a function of the mass. In this combination, 
only physical quantities are compared, that do no longer depend on a particular choice of a 
scale at which one fixes the coupling g, as was the case in Fig. [2j While perturbation theory 
clearly breaks down above a certain ratio m/T > 0.3, RG methods and the 2PI approach 
give consistent results up to twice this value. As noticed previously [14J . the perturbative 
g 2 contribution seems to be a surprisingly good approximation for the behavior at larger 
m/T - but only for the plotted ratio (P/Pq) versus (m/T): the g 2 curve representing P/Pq 
in Fig. [3] is identical to the corresponding curve in Fig. [2] (right panel), since, as already 
mentioned, at leading order in the coupling constant, m/T = g. Thus the reason the curves 
corresponding to the various non-perturbative approximations are brought closer to the g 2 
curve in Fig. [3] can be traced back to the fact that the mass increases much less rapidly than 
g with increasing g, as clearly visible in Fig. |2] (left panel). Note that the dependence of 
the LPA results on the choice of the regulator remains present, although it is less important 
than in Fig. [2j We shall return to this issue shortly. 

We turn now to more technical aspects of the calculations, namely the dependence of the 
results displayed above on the choice of the temperature, or on the choice of the regulator. 
Consider first the dependence on the temperature, which is measured by its ratio T/A to the 
microscopic scale A. As T/A increases and becomes close to 1, our numerical calculations 
loose accuracy for a variety of reasons, the main one being the following: If T is too close 
to A, 2-7rT may become bigger than A and the whole procedure eventually collapses. One 
indeed assumes that the beginning of the flow is not affected by the temperature (so as to 
use the 4-dimensional integration procedures), and this assumes A > 27rT so that there is 
room for a 4-dimensional flow between A and 27rT, with k, w 27rT the point where thermal 
fluctuations start to contribute to the flow. This limitation is illustrated in Fig. HI for the 
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Figure 3: (Color online) Pressure as a function of the mass. The same curves are shown as in 3, but 
for the scale-independent function of pressure as a function of the thermal mass. The perturbative 
results are shown through order g e logg. 

case of the LPA (the phenomenon would be identical in the BMW approximation): the plot 
of the pressure versus the mass is independent of the temperature as long as T/A ^ 1/20, 
while the curve corresponding to T/A = 1/10, deviates slightly from it. 

The right panel of Fig. H] illustrates the dependence of the results on the value of the 
parameter a in the exponential regulator (I2.20p . This is achieved by repeating calculations 
for a set of values of a in the range [1.5, 5]. As one can see, the results are fairly insensitive 
to the value of a, until the mass reaches a value of the order m/T ~ 0.5 -r- 0.6, where a 
sizable dependence starts to be visible. Note that this is the value of the mass where we 
observed possible cut-off effects when the temperature is too large. In the present case, the 
temperature is not to blame. However cut-off effects may show up in the fact that integrals 
are done with A as ultraviolet cut-off. In most cases this is redundant and does not interfere 
with the regulator (whose derivative with respect to k also provides an ultraviolet cut-off). 
However in the mass region m/T « 0.5 -r- 0.6, the two cut-offs interfere, and produce a 
sensitivity of the results to the regulator: indeed a larger a allows larger momenta. 

To investigate this further, we have considered the other regulator mentioned in Sect. [Til 
namely, Eq. f 12 . 2 1 [) . This regulator has the property to cut-off more efficiently the large 
momenta (because of the presence of the g 4 term). And indeed this is what one sees in 
Fig. [5j The left panel indicates that a higher temperature can be reached without affecting 
the results. The right panel shows that the dependence on the regulator parameter a has 
basically disappeared. A further confirmation comes form the study of the dependence of 
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Figure 4: ( Color online) Pressure as a function of the mass, within the LPA and the exponential 
regulator (|2.20[) . Left: comparison of results for various values of the temperature T/A. Right: 
comparison of results for various values of the exponential regulator parameter a = 1.5,2,3,4,5. 
In the right panel, the curves are extracted for T/A = 1/20 or smaller (so that the plots in this 
panel are unaffected by a T/A dependence). 
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Figure 5: (Color online) Pressure as a function of the mass. Comparison of LPA plots at various 
values of the regulator parameter a = 1.5, 3, 4, 5. 

the results on the temperature T switc h at which one changes from the 4-dimensional to the 
3-dimensional flow: essentially no such dependence is observed with the second regulator as 
long as T sw it c h ~ 47r. One should emphasize however that the difference between the two 
regulators is tiny. Thus for instance, the plots in Fig. [6] show the derivative of the pressure as 
a function of k. The oscillatory behavior at the beginning of the flow is generic for regulators 
that cut-off the sum over the Matsubara frequencies, and the particular oscillations exhibited 
in Fig. |6] may be understood from the analytical structure of the regulator analyzed in detail 
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Figure 6: (Color online) Derivative of the flow of the potential difference between thermal and 
vacuum contribution. Comparison of LPA contributions for the exponential regulator (|2.20p (left) 
or ()2.2ip (right) for various values of the parameter a = 1.5,2,3,4,5. The larger the value of a, 
the more oscillations are introduced through the regulator. 

in Appendix B. One sees clearly that the larger a the larger are the oscillations. On the 
other hand, there is hardly any visible difference between the two panels corresponding to 
the two different regulators. Only if one zooms in a lot, can one see that f)2.2ip is slightly 
better in cutting off dP K at larger values of k than (I2.20p . The plots show curves for g\ = 2.8, 
corresponding to g(2irT) = 1.1 and T = 10~ 2 A (that is, the same parameters as the right 
column of Fig. [1]). 

IV. CONCLUSIONS 



In this paper we have applied a powerful approximation scheme to solve the NPRG 
equations and we have calculated various thermodynamical quantities for a scalar field theory 
at finite temperature, in a range of values of the coupling constant covering both weak and 
strong coupling regimes. Of course, because of the presence of the Landau pole, arbitrarily 
large values of the coupling cannot be reached. However the range of values of g that can 
be explored allowed us to demonstrate the stability of the results as successive orders in the 
approximation scheme are taken into account, in sharp contrast to the strict expansion in 
terms of the coupling constant. A similar stability also emerges in other non-perturbative 



schemes, such as "screened perturbation theory" |l9J, or the 2PI effective action 20( (see 
also [2l|). 
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In fact, the results that we have obtained using the BMW approximation differ very little 

n, 

from those that we obtained previously using the LPA |14|]. This may be attributed to 
the fact that the thermal mass provides an infrared cut-off that reduces the contributions 
to the pressure of the long wavelength modes, which are the most strongly coupled. The 
thermal mass also provides a threshold that hinders the mixing with complex multiparticle 
configurations, making the quasiparticles well defined. This latter aspect is explicitly verified 
by the small deviation of the field normalization from unity, as obtained within the BMW 
approximation. These two effects conspire to give the approximation scheme a remarkable 
stability, and contribute to the robustness of the quasiparticle picture. 

We have also compared the results obtained within the BMW scheme with those of the 
simple 2-loop 2PI resummation used in [14J: both methods lead to very similar results in the 
extrapolation to strong coupling. This is not too surprising, given the agreement already 



observed between the LPA and the 2PI method in Ref . 14 1. The stability of the 2PI scheme 
itself can be assessed from the 3-loop calculation of Ref. 20j. The detailed comparison be- 
tween the latter calculation and ours is not straightforward however because of the different 
renormalization schemes used in the two cases. However, the main message is essentially the 
same: the main qualitative difference between the 2-loop and 3-loop calculations is the pres- 
ence of momentum-dependent self-energies in the latter, in contrast to simple mass terms in 
the former. The small difference observed between the results in the two cases, corroborates 
the conclusions that follow from our NPRG analysis about the robustness of the quasipar- 
ticle picture. In fact the stability of the results lead us to conjecture that any corrections to 
them are presumably very small. It would be interesting to have lattice calculations allowing 
us to test this conjecture for values of the couplings where perturbation theory breaks down. 
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Figure 7: (Color online) Integration domains for thermal and vacuum pieces. For the zero- 
temperature piece, Euclidean invariance imposes a four-dimensional sphere for a given cutoff A 
(blue dashed circle) , while the summation over Matsubara frequencies up to a maximum frequency 
independent of the three momentum, implies a four-dimensional cylinder (red solid square). For 
finite A, this mismatch can give rise to spurious contributions to the pressure. 

Appendix A: Numerical integration of the pressure flow 



The calculation of the pressure involves a delicate cancellation between thermal and vac- 
uum contributions at large momenta. The thermal pressure is obtained from the minimum of 
the effective potential, after subtracting the vacuum contribution, as indicated in Eq. ( 12. 8p . 
It is of the form P = J ' dn {p K V^ — d K V^ =0 ), where the vacuum (T = 0) integral with 
vacuum propagator is subtracted from the thermal sum with thermal propagator: 



dtVT - 6W T=0 = 4T 



"max 

E 



d 3 q 1 

——^d t R K {q)G T K {q,p)-- 



d 4 q 



d t R K (q)Gi=»(q,p) 



(Al) 
with f d 3 q = f 2<A 2 t/ 3 q and f d 4 q = J 2 2<A2 d A q. As can be seen in Fig. [7J the two 
integration domains do not match. For large cutoffs A, this mismatch may give rise to 
spurious contributions to the pressure. 

One solution could seem to be to restrict the range of the Matsubara sum to a 4- 
dimensional sphere. This turns out to be problematic for the following reason: Good con- 
vergence properties are only obtained if sufficiently many Matsubara terms are summed up, 
but at large values of q, only few Matsubara terms would contribute if the restriction to a 
4-dimensional sphere were imposed. 

It is better then to fix the mismatch in an alternative way, i.e., by extending the vacuum 
integration domain so that it matches the domain of the Matsubara summation. That is, 
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we write 

" A d 3 q 



1 " max /"A Ji 

dtVj - d t V^ = '-T £ y ^9A(g)Gl(^) 

" J /*""" ^ / ^WGTIW), (A2) 

2 i-TO.max 27r i ( 27r ) 



90, max 

for go.max = 27m max T. Practically, only the pressure is sensitive to this correction. At each 
integration step in t direction, the vacuum integral is calculated twice: Once as in Eq. ( lAip 
to follow the vacuum flow, and additionally as in Eq. (IA2[) to obtain the contribution for 
the pressure. The values of G^ =0 (q,p) that are only known on a grid in (|g|,p) coordinates 
have to be interpolated to obtain values on a three-dimensional grid (go, |q| , p) . Since for 
each temperature the switching point between 4D vacuum integration and 3D Matsubara 
summation varies, the vacuum piece for Eq. (1A2|) has to be calculated for each temperature 
separately. 

Appendix B: Oscillatory behavior of the flow 

As we have seen in Sect. IIII1 the flow of the pressure exhibits an oscillatory behavior at 
the beginning of the flow. This generically occurs whenever the sum over the Matsubara 
frequencies is limited to a finite number of terms. We show in this Appendix that the 
oscillations seen in Fig. [6] at the beginning of the flow can be understood from the analytic 
structure of the particular regulator that we are using, and can be simply calculated for 
small values of k. The analysis is performed within the LPA with the exponential regulator 
(12.201) with a — 1. (Note that much wilder oscillations than the one discussed here are 
observed when one uses a "hard" regulator, such as the Litim regulator |22j.) 

Instead of performing the sum over Matsubara frequencies in Eq. (12.61) explicitly, we can 
convert it into a contour integration in the following way: 

00 1 rioo -1 

T J2 /(*> = *"*) = 7T-. / dq - [/(go) + /(-go)] (Bl) 



too 
ioo+e 



dqo [/(go) + f(-qo)} 



2m J_ ioo+e ™ lJ vw/ J v W/J eVT _ 1 • 

The resulting integrals can then be calculated by closing the vertical contours with semi 
circles of infinite radii, and summing over the encircled residues. 
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Figure 8: Contour plot in the complex qo plane of the integrand of Eq. (|2.6j) . for |q| = 1, m = 0.01, 
and k = 0.5. The left bow of poles corresponds to trivial poles whose residues can be calculated 
analytically. The right bow of poles can only be obtained numerically. 



It is convenient to group the poles of the integrand in Eq. (12.61) into two categories (see 
Fig. [S]). There are the poles associated with the regulator, which we shall call trivial poles, 
and those of the propagator, that is the zeros, of 



G K 1 (q,p) = q 2 + m 2 K + R K {q). 



(B2) 



The positions of the trivial poles are determined by 



exp[(-g 2 + q 2 )/«: 2 ]-l = 0. 



(B3) 



This is satisfied for (— q 2 , + q 2 ) /k 2 = 2irin, with n G N, with the solution q = 
_t\Al 2 — 2ninK 2 . For |q| = 0, these poles lie at (for Rego > 0) go — k(1 ±i) y/rm. For 
|q| 7^ 0, the poles lie along a hyperbola. For completeness, we give the explicit real and 
imaginary parts of go : 



qo 



Tin) k 4 H 

4 



q^ 



±ii 



Tin) k 4 H 

4 



2 



(B4) 



The non-trivial poles can only be found by numerically by solving 



- q 2 Q + q 2 + ml + R K (-q 2 + q 2 ) = 0. 



(B5) 



It turns out that the non-trivial poles lie to the right of the trivial poles in the complex 
plane. Due to the statistical factor in Eq. ( IB II) . the corresponding residues are therefore 
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exponentially suppressed. In fact, the first few trivial poles are sufficient to accurately 
approximate Eq. (12. 6p for large values of k. 

Let us then calculate the residue for the integrand of Eq. ( 12.61) at the trivial pole positions 
(IB4I) . First we note that we can write the derivative of Eq. (I2.19P for the regulator ( 12.201) as 



d t R K ( q ) = -{<i 2 RM + R 2 M- 



(B6) 



Since R K (q) diverges at the trivial pole, the propagator (IB2I) is dominated by the regulator, 
and the residue is independent of the mass m K . We have 

1 -,„_.._. _ 1 2q 2 



Res , ,-^JT T^dtRMG^q.p) = Res — -jf — -^R K {q) 



q 



2 -qi(n) + q 2 



e qo(n)/T _ I - q ^ n j 

with the poles given by qo(n) = ±A/q 2 — 2ninK 2 . The result follows easily from 



(B7) 



Res R K (q) 

qo^qo(n) 



-K 



2 -9oW + q 2 



2q (n) 



(B8) 



The oscillatory behavior follows from the distribution function l/(e qo ' T — 1) with complex 

number q . 

The thermal contribution is directly given by the second line of (IBlj) . One can write the 

contribution of the trivial poles as (instead of adding the residues above and below the real 

axis, we simply can take two times the real value of one of the residues; another factor 2 

comes from /(go) + f(— <?o); the factor 4tt = j dQ) 

1 
dq 2Re — i- 

or explicitly: 



dt [VI {p) - VT*(p)] - 2^-g- £ r dq2Re 

^ ' n>l 



e qo(n)/T 



, - gg(n) + g 2 
qo(n) 



-q 



(B9) 



d t [vl {p) - vI=\ P )] 



1 

2^2 



n>l 



Y^ / dq2Re 



2iririK 



q 



?\/q 2 -2iTiK 2 /T _ j ^Jqi _ 2-KinK 



(BIO) 



For large k this gives a good approximation to the thermal pressure contribution, as can 
be seen in Fig. [9j Taking only the first term n = 1 already gives a good approximation for 
large values of k. 
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Figure 9: (Color online) Numerical solution (LPA) of the thermal contribution of Eq. (|2.6j) . com- 
pared to the first three terms of (JBIOD with n = 1, 2, or 3, as a function of t = Inn, for T = A/20 
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